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1. Introduction 

It is well-known tliat quantum effects lead to the loss of charge by charged black holes 
P, [21 [3], H], and that this phenomenon on one hand is independent on the fact that 
there can be a contextual evaporation phenomenon (i.e. also extremal black holes, with 
zero temperature, are involved in this spontaneous loss of charge), on the other hand its 
being related to the Schwinger instability of vacuum in presence of a constant electric 
field has been pointed out. The latter topic can be brought back to the more general 
class of phenomena of quantum instability in presence of an external field, which still 
has in the Schwinger calculation the most relevant and important contribution [5l [6]. 
See also [3 [8]. 

In previous studies devoted to this topic in the case of black hole backgrounds, we mainly 
focused our attention to black holes of the Kerr-Newman family [9], also in presence of 
a cosmological constant [lOl [TT]. In the latter case, we were able to perform an exact 
calculation for charged Dirac fields in four dimensions in three special cases [12]: the 
ultracold I, ultracold II and Nariai charged black hole backgrounds. The calculations 
were performed both in the so-called transmission coefficient approach and in the zeta- 
function approach, obtaining so a double check for our calculations. 
Herein, we complete our analysis performed in [12] by taking into account the case of 
charged scalar fields, and provide an exact calculation for their instability on the given 
black hole backgrounds. We recall that in similar calculations one is far from being able 
to reach exact results (e.g. in the Reissner-Nordstrom case only a WKB approximation 
is available). We also point out that, contrarily to what one could naively expect, the 
scalar field analysis presents to some extent more difficulties than the analysis for the 
Dirac field, because of some mathematical subtleties occurring in the scalar case: we 
limit ourselves to mention the (open) problem of a rigorous mathematical setting for 
the Klein-Gordon equation minimally coupled with an external electrostatic potential in 
presence of event horizon(s), requiring an analysis involving the so-called Krein spaces 
in place of the more standard Hilbert spaces occurring in the analysis of the Dirac equa- 
tion. Still, we can perform with some ingenuity zeta-function calculations and show that 
the imaginary part of the effective action coincides with the one calculated by means of 
the transmission coefficient approach. Another peculiar behavior emerges in the scalar 
field case when one takes into account the behavior of the field in the ultracold I case: a 
bad behavior at infinity occurs for the wave function, but an analysis in terms of fluxes 
allows to determine the transmission coefficient. Moreover, in the Nariai case, the scalar 
nature of the particle is at the root of the possibility to obtain a change of sign in a 
quantity A (cf. eqn. ( |68ll ) due to the presence of a term — | which is instead missing 
in the analogous quantity for the Dirac case (see [12]). This may cause a change in the 
behavior of the imaginary part of the effective action, as we shall see. 

In the cases ultracold I and Nariai, which are involved with a non-zero background 
temperature, thermal effects on the quantum instability are also considered. 
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The plan of the paper is the following. In section [2] we sum up some aspect of 
the transmission coefficient approach which are relevant for our paper, and then extend 
our analysis [12] concerning instability of thermal state induced by the pair-creation 
effect to scalar fields. In sections [3l [Hand [5] we take into account the cases of ultracold 
II, ultracold I and Nariai charged black hole backgrounds respectively. In section [6] 
conclusions are drawn. 



2. Vacuum instability and Thermal state instability 

We discuss in the following some aspect of the problem of vacuum instability and of 
thermal state instability induced by it. The Dirac case was discussed in jl2] . 



2.1. Vacuum instability 

For completeness, we summarize some aspect of the transmission coefficient approach in 
the case of scalar fields, following [3l [131 E] • We are mainly interested in the probability 
of persistence of the vacuum. Let us introduce, for a diagonal scattering process [3] 

nr = i?.nr^ + T.pr^, (1) 

where rij stays for a negative energy mode and pi for a positive energy one. Tj is the 
transmission coefficient and Ri is the reflection one. Moreover, as in [3], we deflne 

V^ ■■= \Ti\^- (2) 
Then, it is possible to show that for bosons one gets 

\R^\^ = l + V^, (3) 

which accounts for the well-known superradiance phenomenon. 
The persistence of the vacuum is given by [3] 

Po = llP^,o = e-''-^, (4) 

i 

where ptfl is the probability to have zero pair in channel i, and then 

2ImW = J]log(l + m) = Yl (5) 



k 

k=l 



2.2. Thermal state instability 

We have to take into account that, in the case of the ultracold I manifold and also in 
the Nariai charged case, there exists an intrinsic thermality of the background manifold 
which is associated with the presence of non-degenerate horizons. As a consequence, 
the real quantum state to be considered is not the vacuum state in a traditional sense 
(i.e. absence of particles), but the thermal state associated with the aforementioned 
temperature (we recall that we are dealing with special manifolds endowed with a single 
temperature even if two different non-degenerate event horizons are involved). As a 
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consequence, we construct the Hartle-Hawking state for our thermal geometries, by 
adopting the same attitude as in [T2] . 

We point out that the following construction holds true in general, even if we are 
interested in it for our specific analysis. We adopt the thermofield dynamics formalism, 
and define a thermal state \0{f3) > characterized by an inverse temperature /?. This 
state is annihilated by suitable operators ai{P), ai{P), bi{P),bi{P) (and conjugated ones) 
which are labeled by a complete set of quantum numbers / and is related to "standard" 
annihilation-creation operators a/, di, bi, bi (and conjugated ones) via a formally unitary 
transformation: 

ai = c+{(3)ai{P) + s+{di)\(3), (6) 
k=cTm{P)+sT{ky{(3), (7) 

and analogues for hermitian conjugates, with 

4 = ^ (8) 

S;" = ^ =, (9) 

where 4'~^,<f)~ stay for chemical potentials for particles and antiparticles respectively. 
Moreover, it holds 

(ci? - {si? = 1, (10) 

{c^r-{s;r = i. (11) 

We also introduce standard Bogoliubov relations between "in" and "out" operators as 
follows: 

ar = piar+T;{br)\ m 
br = pibr+T;{ary, (13) 

where |pip — |T;p = 1. Compare also [15]. Note that we limited ourselves to consider 



diagonal transformations, as it is the interesting case for our considerations. 

In order to check how thermal effects affect instability of quantum fields, we 
consider, in place of the usual < m|(af"*)''^af"*|0 in >, which gives the number of 
out-particles on the in-vacuum, the following quantity: 

iV+, :=< 0(/5) zn\{aryar'W) m >, (14) 

and check if deviations from pure thermality appears in the distribution. Equivalently, 
as in [151, we can define 



iV+, :=< m ^n\{ar')^a^' - {a'n^o^rW) ^ >, (15) 

which just signals us the deviation part (or it is zero). 
It is easily shown that 

Ntut = \Ti?[{stf + {sif + l\. (16) 



Quantum instability on charged Nariai and ultracold black hole manifolds 



5 



where |T;p is the transmission coefficient. When (p^ = (p = 0, as in the case of our 
black hole background, we obtain 



coth|^'"--^'Vcotl/^'l-'l+^' 



which is easily realized to coincide with the result displayed for the boson case in [T5] 
when = 0, and matches the results in [12] for the Dirac case. Note that (fTTl) can be 
used also for the Reissner- Nordstrom case, where the coefficient |T/p is known only in 
the WKB approximation [Il[2]. 

3. Ultracold II case 

The ultracold II metric is obtained from the Reissner-Nordstrom-de Sitter one in the 
limit of coincidence of the Cauchy horizon, of the black hole event horizon and of the 
cosmological event horizon. See [161 [17]. particular, the metric we are interested in 
is 

ds^ = -de + dy"^ + ^{dO"^ + sin2(^)d02), (18) 

with ?/ G M and t G M. The electromagnetic field strength is F = —\/lidt A dy, and 
we can choose Aq = \/~K.y and Aj = 0, j = 1, 2, 3. It is also useful to define E := ^/X, 
which represents the modulus of the electrostatic field on the given manifold. We note 
that it is uniform, and then one expects naively to retrieve at least some features of 
Schwinger's result, as in the Dirac case [T2]. 

3.1. The transmission coefficient approach 

Let us consider the Klein-Gordon equation in the given manifold 

[-{-zd, + eEyf - dl - 2kVl + /i^]^ = 0, (19) 

where and e are the mass and the charge of the scalar particle. We assume eE > 
for definiteness. In agreement with the possibility to perform variable separation, let us 
set 

^{t,y,n) = e-''^'YiU^)i,{yy, (20) 
then one obtains the following equation for ip: 

^{y) = {^^^-{uj + eEyr)^{y), (21) 

where fj.f = 2M{1 + 1) + /x^. By defining (cf. [3]) 

e = -^{^ + eEy), 
y eE 

. 1 .2 
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u = V2e-'^^, (22) 

one obtains 

^(0 = (A-eXO, (23) 
whose solution is 

^ = Dn{u), (24) 

which is a parabohc cyhnder function. The calculation is completely analogous to the 
one performed in [3j, and as in [3j one can easily show that the transmission coefficient 
satisfies 

2 

|Ti|2 = e""^ = e"^^. (25) 

The latter expression coincides with the WKB approximation for the same coefficient 
pT] (that calculation is for Dirac particles, but it is easy to realize that for scalar particles 
the result is the same, apart for the obvious replacement k"^ ^ l{l + 1)). This means 
that the WKB approximation is actually exact for the given case. We have the exact 
transmission coefficient. As in [3l [2] we can determine the degeneracy factor, and one 
obtains 

^ = ^ E(2^ + 1) + ^~^)' (26) 

1=0 

where S is the spacetime volume of the (t, ?/)-part of the manifold. 
3.2. The (^-function approach 

We can use the ^-function regularization to compute the effective action. The spectral 
zeta function for the Euclidean Klein-Gordon equation is given by 

as) = / z^-'K{z)dz, (27) 



r(3) .0 

with kernel 

K{z) = T:j:e~^-(9r+ieEyy'-d^-2AS7l+fi'']z^ (^28) 

where r stays for the Euclidean time. To compute the trace, note that the operator 
— 2AVq + /x^ commutes with the Klein-Gordon operator so that it contributes with the 
eigenvalues fif = 2M{1 + 1) + ^'^ with degeneration (2/ + 1). Next, noting that the 
operator p := —id^ commutes with A := —{dt + ieEyY — dy we can restrict on the 
eigenspaces having eigenvalue u for p. Thus, A = + cEyY — dy which describe a 
harmonic oscillator with eigenvalues eE{2n + 1). Independence on u shows that such 
eigenvalues are degenerate so that if D is the degeneration we can write 

CO CO oo [eE~^i^]z 

^(^) = E E + l)/^e-[^^(2"+^)+'^^]^ = J2^il + l)D^^^^^. (29) 

1=0 n=0 1=0 
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We can determine the degeneration factor as done in pKj. We then obtain D = ^S, 
where S is the spacetime volume of the {t, ?/)-part of the manifold. 
The Euclidean action is 

Se = -C'(o), 

and one finds 

W = ItuSl, 

where Sl is the Lorentzian action, in our case obtained by — > iE. Explicitly 



as) = j:"£^si2i + i)l 

1=0 ^ 



2eE 



(30) 



where C(s; a) is the Hurwitz zeta function and 7 is a renormalization scale, henceforth 
put equal to 1 (also in the ultracold I and Nariai case). We also put 



1=0 

After Lorentzian continuation we find 



Wi = ^{2l + l)Re 



+«4log(2.eE) + llog2.-C(0;-.|^ + i) 



(31) 



One can notice that: 
,2 



Re 



i-^i- + - log27r - logT - - 
2eE 2 2 ^ ^ ^2 2eE' 



1 , 

-log(l + e-'^^). 



Thus, the final expression for the imaginary part of the Lorentzian action is: 

eES 1 t^f 

Wi = ^i2l + l)-\ogil + e--^), 

which coincides with 



(32) 



4. The ultracold I case 



A second extremal limit of the Nariai background is given by the type I ultracold solution 
when r_ = r+ = Tc. The metric is |l7j 

ds^ = -x^dij^ + dx^ + ^{de^ + sm\e)d<P^), (33) 

with X ^ (0, 00) and ^/^ G M, and the electromagnetic field strength is F = y/Xxdx A dtjj. 
The spacetime presents the structure of a 2D Rindler manifold times a two dimensional 
sphere (with a constant warping factor). One gets Fq^ = -,rQQ = x^^ls = 

— sin(6') cos(^), = cot(6'). We can choose Aq = ^x^ and Aj = 0, j = 1,2,3 as 
potential. This case is a little bit more tricky than the ultracold 11 and we will adopt 
a different strategy to define the transmission and reflection coefficients. However, we 
will again be able to compare this approach with the zeta function method and the two 
results are the same. The expression of the effective action for a scalar field in this 
background fits our previous result for the Dirac case [T2]. 
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The transmission coefficient approach 



In order to compute the wave functions for a scalar field in this background we search 
the solutions of the Klein-Gordon equation, for the variable x, that for this metric is: 



X, ^ X, 



+ l^t 



* = 0, 



(34) 



where to perform variable separation we pose 0(t, x, Q) = e "^'^y;^(Q)^(x). Making 

2 

the change of variable t — 



* = 0. 



Next, we set, as usual, the wave function in the factorized form: 



*(x(t)) = ei^^¥tF(i), 



(35) 



(36) 



so we obtain the following confluent hypergeometric equation or Kummer's equation: 



td^F + {l + iuj + ieEt)dtF - -{/if - ieE)F = 0, 



(37) 



whose general solution is 



Fit) = a$(^ 



l + icu; -leEt) + /3r''^$( ^ 

2ieE' ' ' ^2 2ieE 



icu; 1 



-ieEt), 



where ^{a;c;z) is the usual Kummer function (or 1-st kind confluent hypergeometric 
function). Then, the general solution is 



*(x(i)) = 62 



ieEt 



at-^'^^i- 



2ieE 



■,1 + iu; —ieEt) 



— icu; 1 — icu; —ieEt) 



(38) 



The asymptotic behavior of the wave function can be determined using for \z 

Tic) 



00: 



$(a;c; z) 



z)-"(l + 0(l/^)) + ^e^z-('^-")(l + 0(l/z)).(39) 



r(c — a) 

In the ultracold I background the electric field vanishes in f = and it grows indefinitely 
for t 00. Thus, for large t a charged particle is subjected to an increasing force and it 
is accelerated toward infinity. For this reason the particle, for t 00, does not behave as 
a free particle. The asymptotic behavior of the Kummer function refiects these physical 
considerations: the presence in the asymptotic expansion of the wave functions of terms 
proportional to shows that the behavior of the particle is far to the one of a free 
particle. Thus, one can not define the transmission or reflection coefficients in the usual 
way. However, we can dcflne them using a slightly different strategy. First we compute 
the Klein-Gordon conserved current. Then, in the region t one expects to find a 
fiux of matter coming from t ^ and also a refiected one from large t. Instead, in the 
region t ^ 00, due to the electric field, one expects just the presence of a transmitted 
fiux, the one started from small t region, and no refiected one. These considerations 
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allow, in a quite straightforward way, to define the transmission coefficient as the ratio 
between the transmitted flux at t ~ oo and the incoming flux at i ~ and the reflection 
coefficient as the ratio between the reflected flux and the incoming flux at t ~ 0. Let 
V = the 4- velocity of the static coordinate observer, j the conserved current and dT,^'^ 
the surface element through which we would compute the flux. Then, the associated 
infinitesimal fiux is = ^^tiupcr3''v'^dYl^^^ . As we are interested to the fiux in the x 

direction (x = e^^), dE^*^ = [(5^5^ — 5^5^) /2]d%, where d% is an infinitesimal surface 
element, and then 

Thus, dropping the unessential factor f/2A, we can define the transmission and refiection 
coefficients by looking at the covariant current only. 

The covariant components of the Klein- Gordon conserved current are: = 
— I \^*D^^ — {D^^*)'^]. As we will compute it for the two asymptotic regions t 
and i ~ oo, we need the expansion of the wave functions for small and large t. For i 
we obtain: 

making the change of variable t— \e^^: 
and restoring the time dependence i}): 

with a = q;(|)5'^ and b = Finally we obtain the following expression for 

x-component of the conserved current: 

J, ^ u {\a\' - \b\') . (40) 

For t ^ oo, using the expansion of the Kummer function and making the change of 
variable as before, the asymptotic behavior of ^(x(t)) is: 

fa cie* 4 ^ =r + ' 4 e 
with ki = ^- k2 = f and 



2eE' 



,,2 ,,2 



r ( i - A 1 



^ 2 2eE J 

+ iuj) 

\2 ^ ''2eEj 



a ^{-leE) 2+'2eB *'^2*^2 2eE)+2 
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Restoring the dependence: \1/ ~ Cic*"*^ -x-ioj?/; _|_ ^ge"'^*^ -a:-«tj?/;_ r^j.^^ ^-component 
of the conserved current is: 



(41) 



Thus, the transmission and reflection coefficients are (we are considering the crossing- 
level region that appears for < 0): 



eE 


Cl 


2 


2 uj 


a 


2 



|6p 



As explained before, to avoid particles coming from t ~ oo, we impose the condition 
C2 = and we obtain: 



a- 



r(i + zcu)r(i-^(^-|fe; 



Til-iuj) 



+ i 



2eE J 



2 2 
Cl = a{ieE) 2 *2ei5 2 2ei5 



-ieE)- 
Vil+iu) 



r(l + .a;)r(i-^(o;-4)^_ 



1^ ^^2 ^2e£;/l 

To obtain the transmission coefficient we have to compute |ci 



|2. 



I 12 2 21^ TT-^^- 

ci r = a — e 2eB 



cosh(7r(a; - |^)) - e^- cosh(7r|^; ' ' 



sinh(7rco') cosh(7r(co' — 2^)) 

Finally, for the coefficients \Ti\^ and \Ri\^ we obtain: 
i2 „-7r^ sinh(7rco') 



IT, 



— TT — 

— e 2eB . 



I ft|2 = e-"'^- 



cosh(7r(w - 2^)) 

2 

2eE) 



cosh(7r— 



cosh(7r(w- 1^)) 



Observe that |-R/p — |T;p = 1, as expected for bosons. 

As in the Dirac case, the level- crossing region, assuming eE > 0, is determined by 
uo < Q. Pair production is expected to happen only in this region, thus, for eE > 0, we 
must calculate (cf. [3T|) : 

Wi = lj2^ogil + \Tiiuj)\'), 
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for to < 0. We have ^ J dcu (cf. 

An easy computation shows that: 



), where T stays for a finite time intervaL 



log(l + |T,| 



log(|i?. 



log(l + e"''^) - log(l + 



and we have to evaluate the integrah du log ( 1 + e 



In strict analogy with [12] we obtain: 



eE 



4Lb 



W, = IL 

2 27r 



du log ( 1 + e 



eE 



+ :^Li2 



7T^ 

eE 



(42) 



The factor ^ ^ c/cu j amounts to a degeneracy factor and the same geometric 
considerations done in [12] allow us to evaluate it following |19]. The degeneracy factor 
for the scalar case is the same as in the Dirac case and its value is ^ ( f'' duj] = 
eES/2TT, with S = TL where T and L are the sizes of the space time box over which E 
is non vanishing. This value is exactly the same as the one obtaining in the (^-function 
approach. The final result (H2!) for the imaginary part of the effective action coincides 
with the result (!62l) we will find using the zeta function approach. It is worth mentioning 
that the above background implements the physical model analyzed in [20J, apart for 
the fact that in [20] one deals with a 2D model and a further parameter a appears 
(which in our case is equal to 1). The fact that all our geometries allow a Kaluza- 
Klein reduction (compare the discussion in [12]) explains why a correspondence with 
a 2D model is found: the only substantial difference is represented in our case by the 
presence of an effective mass which is given by fif = fi'^ + 2M{1 + 1) replacing the mass 
/i^ of the aforementioned 2D model. 



4-2. (^-function approach 

Also for this background we analyze pair-production with the zeta function method. 
This technique confirms the results obtained with the transmission coefficients approach. 
The Euclidean Klein-Gordon (KG) operator on ultracold I is: 



KG 



-^dl - -9,(x9,) - 2AV^ + + 2ieEUr + {eEf\x'- (43) 



X X 

In the eigenvalue equation KGcp = Xcj) we put: 

= e-*-^ll„(fi)^(x), (44) 

which leads to variable separation, where Yijn{^) are the usual spherical harmonics 
appearing in every problems with spherical symmetry. Then we obtain: 

2 



A. 



l't + ^[^ + eE-r 



We also introduce t 

\2 



^X^? aiid then we obtain: 



l{>^-l^l)t-\iuJ + eEtf 



ij = o. 



(45) 



(46) 
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and introducing z = eEt, we obtain the confluent hypergeometric equation: 

zd'.g + (1 - a; - z)d,g - ^ (^1 - g = 0. (48) 

We require that solutions ^ belong to -^^^[(0, oo), ^] (the measure is inherited from the 
one of the usual scalar product for scalar particles) . It is easy to reahze that this requires 
to consider different solutions for a; < and for a; > 0. Let us flrst consider: 

m = ^iUl- , 1 - eEt\ . (49) 



.2 V e£; , 
We need the quantization condition: 

with n e N, and then: 

K,i = {2n+l)eE + fil (51) 

We have ^ G L^Q, oo), f ] iff < 0. 
The solution: 

git) = (l - + c^, 1 + c^, eEt^ (52) 

requires a further quantization condition: 

and then 

An,;,.. = (2n + l)eE + 2eEuj + /xf, (54) 

and wc can conclude that ip G i^^[(0, oo), dz] iff a; > 0. 
We obtain that the heat kernel 

K{s)^Y.^2l + l)ki{s) (55) 

receives different contributions from different ranges for uj. In particular, for < there 
is, as in the ultracold II case, a degeneracy in uj to be determined, being A^^; independent 
of UJ in that region. We get 

his) = ^e-'^^^e-^^^-^, (56) 

where formally 

D = I dcu. (57) 

J —oo 
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We determine D as in |18j, by comparing the expansion of ki{s) as s 
heat kernel expansion. We obtain 
cES 



D 



27r 



0+ with the 



(58) 



where 5* is the volume of the first 2D factor of the metric. 
We obtain 

C,W = f^(2e^)-V(i(l + |)..^) 

By rotating eE i-^ ieE and looking for lm('i{0), we obtain a first contribution from the 
a; < region which is easily realized to be the same as in the ultracold II case and a 
further contribution from the a; > region which is given by 



(59) 



|i?eC.(l(l 



ieE" 

as to the first term we get 



l) + (log(2e^)-l)/mCi/(i(l + |^),-l)-/mG(^(l + |^; 



— Li2(— e ) _ Li2(— e ) 



the second one is zero, whereas the third one is 



57r 



2 2 

Li2(— e~"^) — Li2(— e"^) 



As a consequence, we obtain 

T ^^^^ 



log(l + e^~) 



T 1 

2tx Ah 



Li2(-e"~) 



-1); 



(60) 



(61) 



(62) 



2n 2 

which leads to a full accord with (H2l) . 

As to thermal effects, in this case we limit ourselves to point out that eqn. ([1] 
holds, but with a pathological behavior associated with the fact that the chemical 
potential is ill-defined unless a spatial cut-off is introduced at x = Xo < oo. The 
same phenomenon affects Dirac particles [T2] . 



5. Nariai case 



We now consider the more general case, that is the electrically charged Nariai solution. 
The manifold is described by the metric [161 EH llZI 



ds' 



i(- sm\x)d^P' + dx') + ^{dO^ + sm\9)dct>^) 



(63) 



with G M, X e (0, vr), and the constants B = ^ - \J I - 12f^^ . 



A= - B aie 



TT. This manifold 



such that ^ < 1, and Li^ :- 
differs from the ultracold 
version, it corresponds to two spheres characterized by different radii. One finds the 



^ ^ ^, ^.^^ ^ \- '^^^ black hole horizon occurs at x 
differs from the ultracold cases because it has finite spatial section. In the Euclidean 
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■^2 
33 



following non-vanishing Christoffel symbols T^^ = cot(x),rQQ = sin(x) cos(x), Tg 
— sin(6') cos(6'), T^^ = cot{9). For the gauge potential we can choose Ai = —Q^ cos{x)^i- 
Also for this more complex case we study pair-production making use of the transmission 
coefficients and zeta function approach. Again these two different methods give the same 
results. For the Nariai case the zeta function approach requires some mathematical 
techniques recently developed in [22] and their application is strictly analogue to the 
Dirac case, exhaustively analyzed in [T2] . 



5.1. Transmission coefficient approach 

We perform variable separation and set (j){ip,x,^) = e~*^'^Yjm(^)^(x); moreover, we 
define f^f = ^ + ^1(1 + !)• We need to find the solution of the Klein-Gordon equation 
for the variable x'- 



15 1 

■—^{uj - eg— cos x)^ - 9^ (sin x^;^) + /if 

sm X A smx 

Let us first change variable, t = — cos x- Then 

1 . 



^(x) = 0. 



1 - t^W - 2t^' 



LU + eQ—ty - fif 



' A 



^ = 0, 



(64) 



(65) 



where the prime is the derivation w.r.t. t. Note that this equation is invariant under 
{t —>■ — t, Q —>■ —Q} so that we can look at the singularity in t = 1 only and obtain the 
properties of the singularity in t = — 1 by Q ^ —Q. Now, near t = 1 



^ 2(1 -t)^" -2^' + 



which has solution \1' = (1 — t)='=2('^+^'3 a). This suggests to set 
m{t) = (l-t)'+(l + t)'-$(t). 



(66) 
(67) 



so that the equation for the function $ is 

(1 - t2)$" _ 2{t - /+(1 - t) + /_(1 + t))$' - [/i^ - + /+ + L - (/+ - /_)']$ = 0. 
Let us introduce 

We are interested in the level-crossing region, which is, for cE > 0, 

-eE <uj <eE. 

In this region one obtains 



li = -{eE±uj) 



We define also 



A = /.f + {cEf - 



(68) 
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Note that the sign of A is not ensured to be positive. To be precise, we should also 
indicate the dependence of A on /, by writing e.g. A^, but, in order to simplify the 
notation, we leave implicit this dependence. Note also that, if /i^ + {eEY — | < 0, then 
for sufficiently high values of / the quantity A passes from negative to positive values. 
The sign of A is associated with a different behavior of the transmission coefficients and 
then of the imaginary part of the effective action. A little consideration allows to draw 
the conclusion that the behaviors in the two different regions (positive and negative) 
are linked each other by analytic continuation. 

We first consider the case A > 0. The general solution of this equation in the level 
crossing region is easily found to be 

= C+F (ieE + J + i\/A, ieE + \- i^/A; i{eE + a;) + 1; ^ ~ 

(69) 



V 2 ' 2 ' ' ' 2 J 

+ C^F (ieE + ^ + iv^, ieE^]-- iVA;i{eE - cu) + 1; ) 



where Fia^ b; c; z) is the usual hypergeometric function. 
We can use the well known relation 

F{a, b; c; z) - r(c)r(c - a - 6) ^. ^ ^ ^ ^ ^ _ ^. ^ _ ^) 
r(c — a)r(c — b) 

r{c)r{a + b-c) _ 

F{c - a,c - b;l - a - b + c;l - z) (70) 

to look at the asymptotic behavior of near the boundaries. 

Let us now introduce the coordinates x = log tan ^. In this coordinate x ~ 0, tt 
become x ~ —oo, +oo and setting 

r]{x) := ^(x(x)), (71) 

T{l + i{eE-u;))T{-z{u + eE)) 

<^ ■— — ; 7= (72) 

T{^-iuJ-i^/A)T{^-iu + iVA) 

^ Til + tieE-u))m^ + eE)) 

T{l + ieE + iy/A)T{^ + ieE -iy/A)' 

^, T{l + t{uj + eE))T{t{uj-eE)) 

^ r{l + i{u + eEj)Tizi-uj + eEj) 
T{l + ieE + iVA)T{l + ieE -iVaY 

we can write 

ri{x) ^ e*(-+-^)-[C_ + C+a] + e-*("+"^)"C+/? x ^ +oo, (76) 
7]{x) ^ e/("-^-'")^[C+ + C-a'] + e-^("-^-^)^C_/5' x ^ -oo. (77) 

If we are searching for the transmission of a particle coming from x ~ — oo, then 
at t — > +00 we must find only the transmitted particle at a; ~ +oo, with positive 
momentum. With the chosen condition, the positive momentum is {uj + eE) so that we 
must set C+ = and Cout = C*- is the coefficient of the outgoing particle. At a; ~ — oo, 
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the coefficient of the ingoing particle is then Cin = C^a\ so that the transmission 
coefficient is 

' Cin a' T{l + i{uj + eE))T{i{uj-eE)) ' ^ ' 

Then, by using known formulas for the Gamma function, one obtains 
~ 2 eE — u sii\h{Ti{ eE — uj)) smhiiri^ eE + uj)) 

\Ti\ = 1= 1= — . ('9) 

eE + UJ cosh(7r(cj - V A)) cosh(7r(cj + V A)) 

We used a different notation for Ti because it is not yet the transmission coefficient such 
that I-Rip = 1 + l^/P- The latter is obtained by noticing that 

\Ti\^ = --|fl|2, (80) 

q 

where r := u + eE and q := uj — eE. Compare also [23j, where a fine discussion upon 
the topic of the Klein paradox is given. 
As a consequence, we find 

|y^|2 ^ sinh(7r( eE - u)) sinh(7r( eE + u)) ^^^^ 
cosh(7r(ci; — "\/A)) cosh(7r(u; + "\/A)) 
In the limit eE ^ uj one finds 



,,2 1 



which, apart for the term — i, is the result which can be obtained in the WKB 
approximation. 

It is easy to show, in the case A < 0, that the only change consists in the replacement 
^a/A I— > a/I A| in the above formulas for the solution and also for a, a', /3, As a 
consequence, we find the following result: 



|r(- 


\+iuJ-y/\ 






{\+iuJ+ y/\ 






|r(l + i{uj + eE))|2|r(z(cj - eE))|2 



' " |r(l + z(cu + eE))|2|r(z(c^-eE))|2 ' ^ ' 

the denominator is the same as in the case A > 0. As to the numerator one finds 



'2 2 cos(7r2;i) cos(7rz2) 

where zi := a/| A| + iuj and Z2 '■= a/| A| — iuj and standard relations for the Gamma 
function are used. As a consequence we get 

|^^|2 ^ - a; ^ sinh(7r( eE - u)) sinh(7r( eE + a;)) ^^^^ 
eE + u cosh(27ru;) + cos(27rA/|A[) 

We need to calculate 

iy=^5^1og(l + |T,|2); (85) 

i 

we do not perform the sum over /, and then we calculate: 

Wi = (2/ + 1)1 5^ log(l + \Ti{u)\'). (86) 
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Let us start from the case A > 0. We have to perform the following integral: 

"■"^ , , cosh[27reEl - cosh[27ra;] , 

duo log{l H ^ -) 

i-eE cosh[27rv A] + cosh[27ra;] 

where the dependence on I is implicit in A; the integral can be rewritten as follows: 



/ := 



(87) 



2e£; log (cosh [27rv/A] + cosh[27re£;]) - II, 



where 



II : 



eE 



(iti;log(cosh[27rv A] + cosh[27ra;]) 



eE 



2neE 



dy\og{p + cosh[y]), 



1 

which is formally the same integral as in the Dirac case. Then we find 

T 



1 r 

+ — 

47r 



Wi = {21 + 1)2^ (eElog 2 (|cosh[27rV^] + cosh[27re£;]) 

-Li2(-e-''^(^+^^)) + Li2(-e'"(^+^^)) 

Li2(-e2-(^-^^)) + Li2(-e-2-(v^-^^))]) . 

As to the case A < 0, it can be obtained by the replacement \/A 
particular, if //^ + {eEY — ^ < 0, one finds that there exists Ic such that 



Klc 



l>lc 



where 



Wi = ^(2Z + 1) (eElog 2 (|cos[27rvlA]] + cosh[27re£;]) 
+ 1_ [_Li2(-e-''^(^^+^^))-Li2(-e2-(^\^-^^)) 
+ Li2(-e2-(^\^+^^)) + Li2(-e-2-(*^-^^))] ) . 



(88) 



(89) 



(90) 
jA^. In 

(91) 



(92) 



By taking into account that Li2(2;) = lA2{z), it is evident that the latter expression is 
real. 



5.2. Nariai in the C,— function approach 

The Euclidean Klein-Gordon operator for the Nariai solution is given by 

A A 

^idr - ieEcosx? dJsmxd^) - BVl + ^? = KG^^\ (93) 

where r = tip. Let us search for the eigenfunctions of this differential operator. We can 
perform variable separation, as usual. Note that BVq and —idr commute with KG^^^ 
and then one can restrict the study of its eigenvalue equation to the eigenspaces of the 
aforementioned operators, i.e. we can write for its eigenfunctions in these eigenspaces 
/(r, X, Q) = e~'^^'^Yi^jn{^)g{x): where yj^^(Q) are the spherical harmonics. Moreover, 
note that 

/i2 B. 



{-BVt,+^')Yl,rr^ = |^fYl, 



(94) 
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The operator KG'^^'>\^j^rn restricted to the above eigenspaces takes the following form: 



A 



A 



„ iiuj — ieE cos x) 
sin X sm x 



.(95) 



Let us introduce the new variable t = —cosx', then the eigenvalue equation for 
KG^^^ L I m becomes: 



g + {l- t')g" - 2tg' = 0. 



(96) 



To transform this equation into an hypergeometric, we set g{t) = (1 + t)'~(l — ty+ilj{t), 
with 



\ |cu±e£;| . 



(97) 



This choice ensures that the solutions belong into the Hilbert space L^(( — 1, 1)) for all 
values of a; G M. As a consequence, equation (|96|) becomes: 

(1 - t^W{t) + [-2t + 2/+(l -t)- 2/_(l + t)] ij'{t) 
X 



+ 



-uj^ -fi^ + j- 2/+/_ +ll-l+ + l'i-L 



m = 0. 



The general solution of this equation is 



Ht) = 2Fi{a+, a_; 2L + 1; ^) + 2^1(0, b; 21+ + 1; ^—^] 



(9^ 



(99) 



with 



a± = i + /+ + /_ ± y ^ - ^2 - /if + ^ + (/+ + /- )2 + (/+ - L )2. (100) 

Note that this solution has a bad behavior in t = ±1. The only possibility to for it to 
lie in ^^((-1, 1)) is that a_ G -N, that is 



7 - ^ + T + 
AAA 



-n, 



n e N. 



Indeed, the spectrum is discrete with eigenvalues 

1 n 2 



(101) 



(102) 



which are degenerate in the azimuthal quantum number m. 

It follows that the heat kernel for the operator KG^^^ is k{s) = his), where 

A;,(s)=Tre-^^^'"' 



(103) 



Actually the sum over uj is an integral due to the continuity of the —idt spectrum. It is 
convenient to split such integration into two parts that are the interval —eE < u < eE 
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and its complement in R. This is because inside the interval the eigenvalues \n,u),i do 
not depend on cj. In this way, we get 



T 

2^ 



2(2/ +1) j^^ J2' 



-l+c.) +^-(ei?) 



duj 



+ {2l + l){eE)J2e 



sA\(n+^+eEY + !^-{eE) 



2_ 1 
4 



The spectral Riemann ^-function associated to the Klein-Gordon operator with kernel 
ki(t) is then 

Us) = ^ [ f"%{t)dt 



T{s) 

T 

2^ 



(2(2/ + l) 



duo 



+(2/ + 1) {eE) 



[n + \ + uy + i^-{eE) 
1 



2 _ 1 
4 



{n + \ + eEy + ^-{eEy 

An analogous computation of one performed in [121 [22] leads to the following expression 
for the imaginary part of the effective action: 

Wi = —{21 + 1) feElog [cosh[7r(VA - e^)] cosh[7r(VA + eE)] 
27r V L 

1 [_Li2(-e-2"(^+^^)) + Li2(-e'"(^+^^)) 



+2eElog2 + 



47r 



(104) 



which coincides with the one obtained using the transmission coefficient approach. The 
same considerations as for the aforementioned approach in the case A < apply in the 
zeta-function approach. 

5.3. Thermal effects 

We find (for definiteness we choose A > 0) 

rout . sinh(7r( eE - u)) sinh(7r( eE + u)) 



< Nr >i 



cosh(7r(co' — V A)) cosh(7r(c<j + v A)) 
X - (coth[7r(cj — </)+)] + coth[7r(|cj| + v^^)]) 



(105) 

with Lf'^ = e{Ao\Tr—Ao\o) = 2eE = ip~ . We recall that in terms of physical (dimensionful) 
variables, by taking into account that Th = and that Wphys = VAuo, in such a way 

that jSphys^phys = 2tXIjJ. 



6. Conclusions 



Our analysis for the scalar case has confirmed the main features we obtained in [T2]: ex- 
act calculations have been performed, both in the transmission coefficient approach and 
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in the zeta-function approach. The latter is more involved but it also provides us much 
more information with respect to the former, indeed the complete one loop effective 
action (and not only its imaginary part) can be obtained by using the zeta-function, as 
known. Differences with the Dirac case are both of general character (indefinite scalar 
product spaces vs. Hilbert spaces) and in particular characteristics of the cases we an- 
alyzed: in the ultracold I case a bad behavior at x = +00, which does not occur for the 
Dirac case, forced us to refer to fluxes in order to compute the transmission coefficient; 
moreover, in the Nariai case, the coefficient A is not ensured to be positive-definite 
(whereas it is positive definite in the Dirac case), and then one is forced to consider 
both cases. It is worth mentioning that this aspect is not new, because an analogous 
problem occurs in the well-known case of the so-called Sauter potential; nevertheless, a 
discussion of that problem for the Sauter potential is often missing (cf. e.g. [24], where 
the case associated with our A < is considered, and e.g. the results in [23] and in 
[Sj (first paper), for Sauter-like potentials, where the opposite case is given). Thermal 
effects, as in the Dirac case, have been shown to affect the discharge phenomenon, with 
a key-role in the pair-creation phenomenon still to be assigned to the transmission co- 
efficient. 
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